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Abstract 

We present a measurement of the time-dependent CP-violating asymmetries in B° — > K^ir decays 
based on 348 million T(4S) — > BB events collected by the BABAR experiment at the PEP-II 
asymmetric-energy B Factory at SLAC. We measure the direct CP-violating asymmetry C K o n o = 
0.20 ± 0.16 ± 0.03 and the CP-violating asymmetry in the interference between mixing and decay 
S K o v o = 0.33 ± 0.26 ± 0.04 where the first error is statistical and the second systematic. On the 
same sample, we measure the decay branching fraction, obtaining B(B° — > K^tt ) = (10.5 ± 0.7 ± 
0.5) x 10~ 6 . All results presented here are preliminary. 
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The recent measurements of the weak phase (3 in b — > ccs decays from BABAR 1 and Belle [2J, 
have reached the precision of the prediction from fits of the unitarity triangle | 3_ , obtained combining 
the information from CP-conserving quantities to the measurements of other CP-violating (CPV) 
processes. The agreement between the two determinations has shown that the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix [I] correctly describes the source of effects in the Standard 
Model (SM). 

With BABAR and Belle collecting more data, one of the major goals of the two experiments is 
to search for indirect evidence of new physics (NP). One possible strategy consists in comparing 
the established value of to independent determinations of the same quantity, obtained from 
penguin-dominated (in SM) b — ► sqq (q = {d, s}) decays 5 

In the SM, the parameters Cf (describing the direct CPV asymmetry) and Sf (describing 
the CPV asymmetry in the interference between mixing and decay) are expected to be consistent 
with the values from b — * ccs decays (namely C/~ and <5/~ sin 2/3). Small deviations from 
this expectation can be induced by additional CKM suppressed contributions to the amplitude. 
On the other hand, additional radiative loop contributions from NP processes may produce large 
deviations. 

In this letter we present updated measurements of the time-dependent CPV asymmetries and 
branching fraction of the decay B° — > K°tt°. The CKM and color suppression of the tree-level 
b — ► suu transition leads to the expectation that this decay is dominated by a top quark mediated 
b — ► sdd penguin diagram, which carries a weak phase avg(V t bV t * s ). If non-leading contributions are 
ignored, the time-dependent CPV asymmetry is governed by sin 2(3. 

The results presented here are based on 348 million T(4S*) — > BB decays collected with the 
BABAR detector at the PEP-II e + e~ collider, located at the Stanford Linear Accelerator Center. The 
BABAR detector, which is described in provides charged particle tracking through a combination 
of a five-layer double-sided silicon micro-strip detector (SVT) and a 40-layer central drift chamber 
(DCH), both operating in a 1.5 T magnetic field to provide momentum measurements. Charged 
kaon and pion identification is achieved through measurements of particle energy loss (dE/dx) in 
the tracking system and Cherenkov cone angle (9 C ) in a detector of internally reflected Cherenkov 
light (DIRC). A segmented CsI(Tl) electromagnetic calorimeter (EMC) provides photon detection 
and electron identification. Finally, the instrumented flux return (IFR) of the magnet allows 
discrimination between muons and pions. 

We reconstruct K® — ► 7r + 7r - candidates from pairs of oppositely charged tracks. The two- 
track combinations must form a vertex with 7r + 7r - invariant mass within 11.2 MeV/c 2 (3.5<r) of the 
etsablished K® mass JB] and reconstructed proper lifetime greater than five times its uncertainty. 
We form tt° — ► 77 candidates from pairs of photon candidates in the EMC that are isolated from 
any charged tracks, carry a minimum energy of 50 MeV, fall within the mass window 110 < m 77 < 
160 MeV/c 2 , and produce the expected lateral shower shapes. Finally, we construct B° — ► K^n 
candidates by combining K® and ir° candidates in the event. For each B candidate two nearly 
independent kinematic variables are computed. The first one is mg, the invariant mass of the 
reconstructed B meson, B rec . The second one is m m i ss , the invariant mass of the other B, B tag , 
computed from the known beam energy, applying a mass constraint to B rec . For signal decays, 
the two variables peak near the B° mass with a resolution of ~ 5.5 MeV/c 2 (~ 31 MeV/ c 2 ) for 
"imiss (ttib)- Both the m m i ss and tub distributions exhibit a low-side tail from leakage of energy 
deposits out of the EMC. We select candidates within the window 5.11 < m m i ss < 5.31 GeV/c 2 
and 5.1294 < tub < 5.4294 GeV/c 2 , which includes the signal peak and a "sideband" region for 

5 Unless explicitly stated, conjugate decay modes are assumed throughout this paper. 
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background characterization. For the 0.8 % of events with more than one candidate, we select the 
combination with the smallest \ 2 = J2i=n° K°( m i ~ rn 'i) 2 / a m i ' where rrii (m^) is the measured 
(etablished) mass and a mi is the estimated uncertainty on the measured mass of particle i. 

The sample of B° — ► K^tt candidates is dominated by random K^ir combinations from 
e + e~ — > qq (q = {u,d,s,c}) fragmentation. Using large samples of simulated BB events, we 
find that backgrounds from other B meson decays can be neglected. We exploit topological ob- 
servables to discriminate the jet-like e + e~ — * qq events from the more uniformly distributed BB 
events. We compute the value of L2/L0, where Lj = £V |p*|| cos#*| J '. Here, p* is the momentum 
of particle i in the T(AS) rest frame and 9* is the angle between p* and the sphericity axis [Oj of 
the B° candidate, and the sum does not include the decay tree of the reconstructed B. In order to 
reduce the number of background events, we require Lz/Lq < 0.55. We also use the distribution 
of this ratio to discriminate the signal from the residual background. Using a full detector simu- 
lation, we estimate that our selection retains (34.3 ± 1.3)% of the signal events. Here, the error 
includes statistical and systematic contributions. The systematic contribution is dominated by the 
reconstruction of K® and 7r°. 

For each B° — > K®tt° candidate, we examine the remaining tracks and neutral candidates 
in the event to determine if the -B tag meson decayed as a B or a B° (flavor tag). We use a 
neural network (NN) to determine the flavor of the -Btag meson from kinematic and particle iden- 
tification information ^U]. Each event is assigned to one of seven mutually exclusive tagging 
categories, designed to combine flavor tags with similar performance and vertex resolution. We 
parameterize the performance of this algorithm in a data sample (-Bfl av ) of fully reconstructed 
B° D^-7r + /p + /af decays. The average effective tagging efficiency obtained from this sample 
is Q = J]c e 5(l ~~ 2w c ) 2 = (30.4 ± 0.3)%, where e c s and w c are the efficiencies and mistag proba- 
bilities, respectively, for events tagged in category c. We take into account differences in tagging 
efficiency (for signal and background) and mistag (only for signal) for B° and B° events, in order 
to exclude any source of fake CPV effects. For the background, the fraction of events (e^) and the 
asymmetry in the rate of B° versus B° tags in each tagging category are extracted from the fit to 
the data. 

Time-dependent CPV asymmetries are determined by reconstructing the distribution of the 
difference of the proper decay times, At = tcp — itagj where the top refers to the signal B° and itag 
to the .Btag- At the T(4S) resonance, the At distribution follows 

e -|At|/r 

V% (At) = x [ 1 ± (S f sin (Am d At) - C f cos (Am d At) ) ] , (1) 

where the upper (lower) sign corresponds to -Bt ag decaying as B° (B°), r is the B° lifetime averaged 
over the two mass eigenstates, Am^ is the mixing frequency, Cf is the magnitude of direct CP 
violation in the decay to final state /, and Sf is the magnitude of CP violation in the interference 
between mixing and decay. For the case of pure penguin dominance, we expect S^-o^o = sin 2/3, 
and C K° s nO = 0. 

We compute the proper time difference At from the known boost of the e + e~ system and 
the measured Az = zqp — -Ztag> the difference of the reconstructed decay vertex positions of the 
B° — > Kg-K and -B tag candidate along the boost direction (z). A description of the inclusive 
reconstruction of the -B tag vertex is given in For the B° — ► K^tt decay, where no charged 

particles are present at the decay vertex, we identify the vertex of the fully reconstructed B using 
the single trajectory from the 7r + 7r _ momenta and the knowledge of the average interaction 
point (IP), which is determined on a run- by-run basis from the spatial distribution of vertices from 
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two-track events. We compute At and its uncertainty from a geometric fit to the T(4S) — > B®B® 
system that takes this IP constraint into account. We further improve the sensitivity to At by 
constraining the sum of the two B decay times (top + *tag) to be equal to 2 r with an uncertainty 

t b o, which effectively constrains the two vertices to be near the T(45) line of flight. We have 
verified in a Monte Carlo simulation that this procedure provides an unbiased estimate of At. 

The per-event estimate of the uncertainty on At reflects the strong dependence of the At 
resolution on the Kg flight direction and on the number of SVT layers traversed by the Kg decay 
daughters. In about 60 % of the events, both pion tracks are reconstructed from at least 4 SVT hits, 
leading to sufficient resolution for the time-dependent measurement. The average At resolution in 
these events is about 1.0 ps. For events which fail this criterion or for which cxa* > 2.5 ps or 
At > 20 ps, the At information is not used. However, since Cf can also be extracted from flavor 
tagging information alone, these events still contribute to the measurement of Cf and the signal 
yield. 

We obtain the probability density function (PDF) for the time-dependence of signal decays 
from the convolution of Eq. ^with a resolution function 1Z(5t = At — Att rue , 0"At)> where At true is 
the true value of At. The resolution function is parameterized as the sum of a 'core' and a 'tail' 
Gaussian, each with a width and mean proportional to the reconstructed OAt, and a third Gaussian 
centered at zero with a fixed width of 8 ps . We have verified in simulation that the parameters 
of 1Z(5t,aAt) f° r B° — > K®tt° decays are similar to those obtained from the -Bfl av sample, even 
though the distributions of OAt differ considerably. Therefore, we extract these parameters from a 
fit to the .Bflav sample. We find that the At distribution of background candidates is well described 
by a 5 function convolved with a resolution function with the same functional form as used for 
signal events. The parameters of the background function are determined together with the CPV 
parameters and the signal yield. 

We extract the CPV parameters from an extended unbinned maximum-likelihood (ML) fit to 
kinematic, event shape, flavor tag, and time structure variables. We construct the likelihood from 
the product of one-dimensional PDFs, since all the linear correlations are negligible. The systematic 
from residual correlations is taken into account, as explained below. 

The PDFs for signal events are parameterized from either a largest sample of fully-reconstructed 
B decays in data or from simulated events. For background PDFs, we select the functional form 
from data in the sideband regions, included in the fitted sample, of the other observables where 
backgrounds dominate 

The likelihood function is defined as: 

-(Ns+N B ) 

C(S f ,Cf,N s , N B ,fs, f B ,a) = x JJ [Nsfst&Pstfuffi] S f ,C f ) + N B f B e c B V B {xi,mS)\ x 

^ S B '' iSw/A* 

II IW ~ fsVsV'six,- Cf) + N B (1 - f B )e B V' B (x t ; a)} , 

iSw/o At 

where fs is the fraction of events with At information (w/ At) and is the fraction of events 
without it (w/oAt). 

The probabilities Vs and Vb are products of PDFs for signal (S) and background (B) hypotheses 
evaluated for the measurements x*i = {ms, m m i ss , L2/L0, cos 9 B , tag, tagging category} and y-i = 
{At,crAt}. In the formula, a represents the set of parameters that define the shape of the PDFs. 
Along with the CPV asymmetries Sf and Cf, the fit extracts the yields N$ and Nb, the fraction 
of events with At information fs and fs, and the parameters a which describe the background 
PDFs. 



10 



200 I — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — f 



CZ3 

a 150 

> 
0) 



100 - 



0) 

W) 50 

* » 



preliminary 



J i i i I i i i I i i i I i i i I i i i I i i i I i i i L 




J I I I I I L= 



5.12 5.14 5.16 5.18 5.2 5.22 5.24 5.26 5.28 

m miS s [GeV] 



Figure 1: m m i ss distribution for signal events on data (dots), obtained using the sPlot technique 12 
to subtract background events. The solid curve represents the shape of signal PDF, as obtained 
from the fit. 



Fitting the data sample of 17058 B° —* Kg-ir candidates, we find N$ = 425 ± 28 signal decays 
with S K o n o = 0.33 ± 0.26 ± 0.04 and C K o n o = 0.20 ± 0.16 ± 0.03, where the uncertainties are sta- 
tistical and systematic respectively. Taking into account the selection efficiency and the number of 
BB pairs included in the fitted data sample, we also obtain B(K°ir°) = (10.5 ± 0.7 ± 0.5) x 10~ 6 . 

Figure n shows the m m i ss distributions for signal events, where background is subtracted using 
the sPlot technique |12j . FigureEJshows distributions of At for B°- and -B°-tagged events, and the 
asymmetry ^4^-o 7r o(At) = [N B o — iVgo] / [Nro + Ng ] as a function of At, also obtained with the 
sPlot event weighting technique. Ng (ATbo) represents the number of events tagged as B° (B°). 

In order to investigate possible biases introduced in the CPV measurements by the IP-constrained 
vertexing technique, we examine B° — > J/ipK® decays in data, where J/tp — ► or J/tp — ► e + e _ . 

In these events we determine At in two ways: by fully reconstructing the B° decay vertex using 
the trajectories of charged daughters of the J/tp and the K® mesons, or by neglecting the J/tp 
contribution to the decay vertex and using the IP constraint and the K® trajectory only. This 
study shows that within statistical uncertainties, the IP-constrained At measurement is unbiased 
with respect to the standard technique and that the obtained values of Sj^ K o and Cj^ K o are 
consistent. 

To compute the systematic error associated with the signal yield and CPV parameters, each 
of the input parameters to the likelihood fit is shifted by ±1<t from its nominal value and the fit 
is repeated. Here, ±lcr is the associated error, as obtained from the -Bfl av sample (for At and 
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Figure 2: Distributions of At for events weighted with the sPlot technique for -Btag tagged as (a) 
B° or (b) B°, and (c) the asymmetry A(At). The points are weighted data and the curves are the 
PDF projections. 
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tagging) or from Monte Carlo. This contribution to the systematic error takes into account the 
limited statistics we used to parameterize the shape of the likelihood. We obtain a systematic error 
of 0.72 events on the yield, and of 0.006 (0.010) on S K o T o (C K o n o). As an additional systematic 
error associated with the shape of the PDF, we also quote the largest deviation observed when 
the individual signal PDFs are floated in the fit. This gives a systematic error of 11 events on the 
yield, and of 0.007 (0.021) on S K o n o {C K o n o). The output values of the PDF parameters are also 
used to associate a systematic error to the selection cuts on the likelihood variables. We evaluate 
the systematic error coming from the neglected correlations among fit variables using a set of toy 
Monte Carlo experiments, in which we embed signal events from full detector simulations. We use 
the average shift in yield (2.3 events) and CPV parameters (0.003 on S K o n o and 0.015 on C K o n o) as 
the associated uncertainty. We estimate the background from other B decays to be negligible in the 
nominal fit. We take into account a systematic error induced on signal yield and CPV parameters 
by this neglected component, embedding B background events in the dataset and evaluating the 
average shift in the fit result: 4.5 events on the signal yield, 0.003 on S K o 7T o and 0.002 on C K o n o- 

For CPV parameters, we evaluate the additional systematic uncertainty related to the fit method 
using the largest difference between the fitted and generated values of S K o n o (0.027) and C^o^o 
(0.003). To quantify possible additional systematic effects, we examine large samples of simulated 
B° — > KgTT and B° — ► J/if)K® decays. We employ the difference in resolution function parameters 
extracted from these samples to evaluate uncertainties due to the use of the resolution function 1Z 
extracted from the -Bfl av sample. We assign a systematic uncertainty of 0.01 on S K o n o and 0.02 on 
C K o n o due to the uncertainty in 1Z. We include a systematic uncertainty of 0.002 on S K o n o and 
0.001 on C K o T o to account for a possible misalignment of the SVT. We consider large variations 
of the IP position and resolution, which produce a systematic uncertainty of 0.004 on S K o n o and 
0.001 on C K o n a. Additional contributions come from the error on the known B° lifetime (0.0022 
on both S K o n o and C K o n o), the value of Arrid (0.0017 on both S K o T o and C K o T o), and the effect 
of interference on the tag side (0.0014 on S K o n o and 0.014 on C K o n o). 

For the branching fraction, systematic errors come from the knowledge of selection efficiency, 
(34.3 ± 1.3)%, the counting of BB pairs in the data sample, (347.5 ± 1.9) x 10 6 BB pairs, and the 
branching fractions of the B decay chain (B(Kg — ► 7r + 7T ") = 0.6895 ± 0.0014 and B(ir° -> 77) = 
0.9880 ±0.0003). [S] 

In summary, we have performed a measurement of the time-dependent CPV asymmetries of 
B° — > KgTT and the branching fraction of B° — > K°ir°. We measured the parameters of CPV 
asymmetry C^o^o = 0.20 ± 0.16 ± 0.03 and S K o n o = 0.33 ± 0.26 ± 0.04, and the branching fraction 
B(B° — ► K°ir°) = (10.5 ± 0.7 ± 0.5) x 10~ 6 . The first error is statistical and the second systematic. 
All the results presented here are preliminary. 
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